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ABSTRACT 

We study the evolution of the mass function in young and dense star clusters by 
means of direct A-body simulations. Our main aim is to explain the recent observations 
of the relatively flat mass function observed near the centre of the Arches star cluster. 
In this region, the power law index of the mass function for stars more massive than 
about 5-6 Mq is larger than the Salpeter value by about unity; whereas further out, 
and for the lower mass stars, the mass function resembles the Salpeter distribution. We 
show that the peculiarities in the Arches mass function can be explained satisfactorily 
without primordial mass segregation. We draw two conclusions from our simulations: 
1) The Arches initial mass function is consistent with a Salpeter slope down to ~ 
1 Mq, 2) The cluster is about half way towards core collapse. The cores of other star 
clusters with characteristics similar to those of the Arches are expected to show similar 
flattening in the mass functions for the high mass ( £ 5M Q ) stars. 
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1 INTRODUCTION 

The mass function of a star cluster changes because of 
both stellar evolution and stellar dynamics. Stellar evolu- 
tion causes the turn-off mass to decrease as the most mas- 
sive stars evolve away from the main sequence, ascend the 
giant branch to ultimately shed their envelopes to turn into 
compact objects. Stellar evolution therefore has a character- 
istic effect on the mass function by truncating it at the high 
mass end. 

The dynamical evolution of a cluster has a more compli- 
cated effect on changes in the mass function. The dominant 
effect here is dynamical friction, which causes the most mas- 
sive stars to sink to the cluster centre on a time scale that is 
inversely proportional to the stellar mass, i.e. the most mas- 
sive stars tend to sink more quickly than relatively lighter 
stars. At the same time, stars less massive than the average 
mass tend to leave the inner regions. As a result of this mass 
segregation, the local stellar population becomes a function 
of the distance to the cluster centre. 
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Mass segregation, though mostly noticeable in the clus- 
ter's central regions, is a global phenomenon. A star cluster 
that is born with the same mass function across its radial 
coordinate will gradually grow a top-heavy mass function in 
its centre and a top-depleted mass function in its outskirts. 
Near the half mass radius, t he mass function remains c losest 
to the initial mass function (|Vesperini fc Heggie|ll997h . 



In this letter, we concentrate on the evolution of the 
stellar mass function in the inner part of young and dense 
star clusters, using A^-body simulations. Our interest in this 
topic was in itiated by the rec e nt accurate meas urements 
published bv lStolte etaLI (|2005l h lKim et all (|2006h in which 
the mass function in the inner ~ 10" from the centre of 
the Arches star cluster was studied. These observations, es- 
pecially the latter, revealed that the mass function of near 
the centre of Arches cluster is a broken power law, with the 
turning point m p ~ 5— 6Mq. We were able to reproduce this 
feature without invoking any special mechanism. Our sim- 
ulations allow us to draw conclusions on the history of the 
dynamical evolution of the Arches cluster. 
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2 DYNAMICAL EVOLUTION OF THE MASS 
FUNCTION 

2.1 Parameters for the simulations 

As a cluster evolves, stars more massive than the mean mass 
(m) tend to sink to the cluster centre whereas lighter stars 
move outwards. For the most massive stars, the time scale 
for dynamical friction is proportional to two-body relaxation 
time, t T : 



(m) 

tdf OC tr 



(1) 



were m+ is the mass of the massive star, which segregates 
inwards. The value of the relax ation time at the cluster's 
half-mass radius, is given by (|Spitzerlll987l . eq. 2.63) 
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Here G is Newton's constant of gravity, M and N are the to- 
tal mass and the number of stars in the cluster and In A is the 
Coulomb logarithm, fo r which we adopt In A = ln(O.OlA) 
l|Giersz fc Heggiell 19961 ). For the central relaxation time, we 
use 



tr, 
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where a so, n and (jn c ) are the three-dimensional velocity 
dispersion, nu mber density and average stellar mass at the 
cluster centre |Spitzerlll987l . eq. 3.37). 

We follow the dynamical evolution of our mod- 
els by means of direct A-body simulations, which 
we carry out with the st arlab software environment 
l|Portegies Zwart et al.l 120011 ). The calculations are per- 
formed on the GRAPE-6 spe cial purpose computer 
jMakino et al.lll997l ; lMakindl200ll l. 

Our numerical experiments are performed with TV" = 
12288 and 24576 stars. For each N, we perform simulations 
starti ng wi t h a f ull range of density profiles for which we 
chose iKinel l|l966h models with the dimensionless parameter 
Wo ranging from 3 to 12. The mass function in our simu- 
lations is described by a power-law, dN/dm = m x , where 
we adopt the Salpeter value for the index (x — —2.35), with 
masses ranging from lMg to 100 Mq. To validate our re- 
sults, we carried out additional simulations with A = 49152 
as well as with a Salpeter mass function with 0.1 Mq as 
the lower limit. It will turn out that the presence of a tidal 
field has little effect on the results, but reducing the lower 
limit to the initial mass function to 0.1 Mq has a profound 
effect on the results, as we discuss below. For clarity we 
mainly focus on the models with 12288 and 24576 stars. 
With these parameters, the relaxation time at the virial ra- 
dius for the 12k models is about 360 A-body time units, 
whe reas for the 24k models t his is 625 A-body time units 
(see lHeggie fc Mathieul l|l98rj fl). 

The c lose proximity of the Arc hes cluster to the Galac- 
tic centre (|Cotera et al.lll992l . [19961 would seemingly require 
the simulations to include tidal effects. And for understand- 
ing the dynamics in the cluster outskirts or the evaporation 



1 For the definition of an N-body unit, or the summary at 
http : //en. wikipedia. org/wiki/Natural_units 



time scale the tidal field will prove crucial. For studying the 
evolution of the central region on the short time scale re- 
ported here, however, the tidal field has negligible effect. 
We support this statement by carrying out additional sim- 
ulations which include the tidal field, and those show no 
discernible effect. We therefore focus on the results of sim- 
ulations without a tidal field. This has the attractive side 
effect that it allows us to scale our results with respect to 
N. We also ignore the effects of stellar evolution. This ap- 
proximation is rectified as on the short lifetime of the cluster 
(2 ± 1 Myr) even the most massive stars remain on the main 
sequence, though some effect of the stellar mass loss at the 
top end of the mass function can be expected. For example, 
a 60 Mq zero-age main sequence star with solar metalicity 
loses about 3 Mq in its first ~ 2.4 Myr jLeieune fc Schaererl 
2001). which has a negligible effect on the slope of the mass 
function. 

2.2 Dynamical evolution towards core collapse 

In our simulations we identify the moment of core collapse 
as soon as a persistent binary forms with a binding energy of 
at least 100 kT (where the energy scale kT is defined by the 
condition that the total stellar kinetic energy of the system, 
excluding internal binary motion, is |AkT). For a cluster 
with a mass function that is consistent with the observed 
mass function in young star clusters, core collapse occurs at 
a more or less constant fracti on of the initial central relax- 
ation time tec ~ 0-2 ± 0. 1t rc (|Portegies Zwart fc McMillanl 
120021 ; iGiirkan . et alj|200lh . In Fig.Q]we plot the moment of 
core collapse as a function of the initial concentration of the 
cluster. The slight dependence of t cc /t IC on Wo, as well as 
the of fset between our results with those of IGiirkan et al.l 
l|2004h is presumably mainly caused by their broader range 
of stellar masses (0.2 < tti/Mq < 120) in the initial mass 
function, whereas here we adopt an initial mass function 
with 1 < m/Mg < 100. An additional effect is expected 
fr om the difference in the number of stars. The simulations 
of IGiirkan et all (|2004l ) were carried out with 10 6 . The sys- 
tematic difference between our results from simulations with 
12k and 24k stars shows that this also affects the results sys- 
tematically. At this moment, however, we cannot quantify 
this effect. 

Dynamical friction causes the massive stars to segre- 
gate to the cluster centre making the mass function flat- 
ter at the higher end, in this region, until the formation 
of a hard binary. In Figure[2]we illustrate the evolution of 
the mass function between r coro and 2r core . We show the 
mass function at birth (top curve, the Salpeter mass func- 
tion), halfway to core collapse and at core collapse (bottom 
curve). We denote the point around which the slope of the 
mass function changes by m p , and the power-law indices in 
higher and lower ends by x m < m and £ m >m p , respectively. 
The apparent decrease in the number of stars in the mass 
function presented in Figure[2] is the result of the cluster 
becoming more concentrated which causes the adopted an- 
nulus (r corc < r < 2r corc ) to become narrower. 

The effect of flattening of the mass function is less pro- 
nounced further away from the cluster centre. This is illus- 
trated in Figure^] where we present the evolution of a; m >m p 
for the 24k simulations with Wo = 5 and for r = to r COTC 
(top curve), for r = r coro to 2r CO rc and for r = 2r corc to 3r core 
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Figure 1. The moment of core collapse (t cc ) in units of the core- 
relaxation time (tic) as a function of Wo- The dashed curve gives 
the results of our simulations with 12k stars and the solid curve 
with 24k. Since we performed only one simulation per set of initial 
conditions no error bars are presen t ed. T he bottom (dotted) line 
denote the results of lGiirkan et alj l|2004h . 
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Figure 2. The mass function of one of the simulations (N = 24k, 
Wo = 5) for an annulus 1 < r/r coro < 2, around the cluster 
centre. From top to bottom, the curves are at times t = 0, t ~ 
0.61f C c, and t ~ 1.05i cc - The solid lines are least squares fits to 
the mass function with a broken power-law. 



(bottom curve). The values of m p , a; m <m p and a; m >m p are 
obtained by a three-point least squares fit to the mass func- 
tion in a predetermined annulus of the simulated date. Note 
that we relaxed the fitting procedure in the sense that the 
mass function is not required to be continuous. The point of 
stalling of the evolution of the mass function can be iden- 
tified by the moment of core collapse, regardless of the ini- 
tial concentration or the number of stars in the simulation. 
Therefore, we normalize the time axis in FigureO to that 
instant. 



Figure 3. The evolution of x m > mp for various radial bins in the 
24k simulation with a Wo = 5 King model. The time is given in 
units of core collapse time f cc . The radial bins are < r / r COIC < 1 
(upper solid curve), 1 < r/r COTC < 2 (dotted curve), and 2 < 
r/r COTe < 3 (lower solid curve). To guide the eye, we plotted 
straight lines through the simulation data. 
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Figure 4. The value of i m > n as a function of distance form 
the cluster centre after core collapse. The solid curve gives the 
average value of x m > mp over the various simulations with 24k 
stars for Wo = 3, 4 up to Wo = 12, the dashed line gives the 
data for the simulations with N = 12k. The squares, bullets 
and triangles with error bars give the obser ved values t aken from 
iFieer et all Jl99Sh . lStolte et al.1 1120051 1 and lKim etafl (I2006TI , re- 
spectively (see Tab.[TJ. The lower thin dashed line gives the value 
of £m>m p for simulations with 24k particles with Wo — 9 and a 
lower limit to the initial mass function of 0.1 Mq. 



2.3 Post-collapse mass function 

After the formation of a hard binary, the mass function 
achieves a quasi steady state. The slope of the high-mass 
end of the mass function varies throughout the cluster. In 
in Figure|4] we show how the mass function for stars with 
m > m v after the moment of core collapse is a function of 
the distance to the cluster centre, being flatter closer in and 
resembling the initial mass function further out. Overplotted 
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Figure 5. The value of m p after the point of core collapse, as 
a function of distance to the cluster centre for various of simula- 
tions. The circles connected with a thin solid line give the results 
for the simulations with 24k stars and with Wo = 5 averaged 
between the moment of core collapse and twice the core collapse 
time. The error bars indicate the variation of the value of m p 
over this time period. More concentrated initial models tend to 
have a slightly lower value of m p , which we illustrate by plotting 
the Wo = 9 simulation as the thick lower dotted line. The upper 
dotted line gives the results of the simulation with N = 12k and 
for Wo = 5. The thin dashed line at the bottom gives the value of 
m p for simulations with 24k particles with Wo = 9 and a lower 
limit to the initial mass function of 0.1 Mq. 

are the observed values of the mass function exponent (see 
Tab.[T] see §[3j) . The results of our simulations with a tidal 
field are statistically identical to those with a tidal field. The 
simulations with a minimum mass to the initial mass func- 
tion of m m j n = 0.1 Mq is plotted as the thin dashed curve in 
Fig. [4] For clarity we did not plot error bars for this figure, 
but the results with a lower limit of 0.1 Mq are inconsistent 
with the observed values. 

In FigureJS] we show the value of m p as a function of 
distance from the cluster centre for various simulations, past 
the moment of core collapse. It turns out that more concen- 
trated initial models tend to result in a slightly smaller value 
of m p whereas simulations with a smaller number of stars 
give rise to a higher value of m p . The behaviour of a: m >m p is 
rather insensitive to the initial concentration of the cluster. 

It may be noted that the results of simulations with 
0.1 Mgas the lower limit of the mass function are not con- 
sistent with the observed values of :r m>mp and m p . 



3 MASS FUNCTION OF THE ARCHES 
CLUSTER 

At a projected distance of about 25 pc from Sgr A*, the 
Arches cluster (a = 17 ,t 45" 5 s , 6 = -27 °49'28" in J2000), 
discovered by ICotera et alj ll 19921 . Il996l ). is peculiar. The 
total cluster mass is about 2 ■ 10 4 M Q . The core radius of 
the cluster (defined as the radial distance from the cluster 
centre where the luminosity pro file drops by a factor two) 
is r CO rc= 5".0 (|Stolte et al.ll2005T ). and corresponds to about 
0.2 pc if we assume that the distance to the Galactic centre 



Table 1. Parameters for the observed mass function of the 
Arches cluster. The first two columns give the range over which 
the mass function is measured, in units of the cluster's core radius 
(r coro ~ 0.20 pc). The third and fourth columns give the range in 
masses for which the exponent of the mass function (last column) 
is fitted. Column 5 gives the ref erence for the mas s fu nction expo- 
nents l: [Stolte et alj J2005h . 2: iKim et all fcOOrJ) , 3: iFiger et al] 
ll 19991) . and the last column the measured value of x between 
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Figure 6. The present day mass function of the Arches star 
cluster between 1 r CO r e and 1.8r CO re. The error bars are taken 
directly from Fig. 5 of lKim et all d2006l) . The three thin lines are 
from Fig.[2]at zero age (top line), at t ~ 0.61t C c (middle line) and 
at t ~ 1.05t cc (bottom line). The thick dashed line is identical to 
the thin t ~ 0.61t cc but then renormalized with -0.35 dex. This 
line produces a satisfactory fit to the observed mass function. 

is 8 k pc. The age of the cluster is 2 ± 1 Myr l|Figer et al.l 
119991) . 

Recently. IKim et alj (|2006h observed the Arches cluster 
using Keck/NIRC2 laser guide star adaptive optics. Their 
observations covered the inner parts of the cluster and some 
control fields at a distance of about 2.4 pc (60") from the 
cluster centre. They subsequently constructed the luminos- 
ity and mass functions down to about 1.3 Mq, in an annulus 
of 5" (about TO r corc ) to 9" (about 1.8 r corc ) from the cluster 
centre. In Table[T]we give the various measurements of the 
mass function and the distance from the cluster centre in 
terms of the observed core radius. 

The data show that the slope in the mass function for 
stars more massive than ~ 5 Mq flattens towards the cluster 
centre (see fig. [4}. For lower mass stars, as well as for further 
out than ~ 4 r corc it is closer to to the Salpeter mass func- 
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tion. The obse rved mass function in the Arches is not a sim- 
ple power-law (jStolte et al.ll2005l ; lKim et alj2006h . We argue 
that in the inner parts of the cluster, r 4r CO re, the mass 
function is best described by two power laws with the break 
around m p — 5-6 M©. The mass function below this break 
(m p ) resembles the initial mass function (x m <m„ = iimf), 
and above m p it becomes flatter (2; m >m p > iimf)- Further 
out than r ~ 4r core the break disappears and the mass 
function becomes gradually better represented with the ini- 
tial mass function. 



4 DISCUSSION AND CONCLUSIONS 

We performed detailed simulations of the evolution of young 
and dense star clusters using direct iV-body simulations, in 
order to constrain the observed mass function within about 
one parsec from the centre of the Arches cluster. The initial 
conditions of our simulations range over the full spectrum 
of King model density profiles. 

The mass function in the central region of the Arches 
cluster is peculiar as it appears to be split in two power- 
laws, one for the stars less massive than 5-6 M© and a much 
shallower slope for the more massive stars. The two power- 
laws fit the observed data between 5" and 9" is marginally 
better (x 2 — 0.5) than a single power-law (x 2 — 0.93). 

The simulations we perform to mimic the Arches cluster 
are able to reproduce this observed broken power-law mass 
function at the observed projected distance from the cluster 
centre (r = 0-4r corc ). The best comparison between obser- 
vations and simulations is obtained if the cluster is about 
half way core collapse (t — 0.4-0.6 t cc ). 

Our simulations, however, are performed without stel- 
lar evolution and without including the effects of an external 
tidal potential. As a result, they are scale-free, and no spe- 
cific choices for the scalings to mass, size and therefore to 
time are obliged. However, the scale-free aspect of our sim- 
ulations hinders the direct comparison to some extend as 
the size scale (in parsec) and time scale (in Myr) are impor- 
tant for an unbiased comparison with the observed Arches 
cluster. 

In the comparison with the observations we adopt the 
same definition of the core radius by projecting the clus- 
ter and assigning luminosities of all stars in our simulations 
using zero-age main-sequence luminositiefl We ignore here 
the fact that very massive stars may become brighter in the 
2±1 Myr lifetime of the cluster, but this only affects the most 
massive stars, whereas the measurements are dominated by 
stars in the mid-range of masses. 

The Arches cluster does not show any evidence for pri- 
mordial mass segregation as our simulations (which were 
initialized without primordial mass segregation) are able to 
satisfactorily reproduce the observed mass function over the 



2 For comparison, we adopt the observers' definition of core ra- 
dius: the point where the surface brightness drops to half its 
central value. This is simil ar to but dif ferent from d efinitions 
used in theoretical works of llSpitzerlll987l. eq. 1-34) and dAarsethl 
2001 eg. 15-4). In sta rlab we adopt the method as discussed by 
Heggie. Trenti. fc Hull 1120061 1. The latter definition of the core 
radius systematically is about twice the observers' definition. 



entire range of observed masses and distances from the clus- 
ter centre. Note also that the presence of primordial gas 
which failed to form stars does not seem to have affected 
the early cluster evolution, as the observed cluster structure 
at an age of 2 ± 1 Myr is satisfactory explained with the 
simulations, which do not include gas dynamics. The initial 
mass function of the Arches cluster is then consistent with a 
Salpeter slope between 1 M© and 100 M© without the need 
for a radial dependence. There seems to be no need for a 
large population of stars less massive than ~ 1 M© . 

In fig.|3] we show the evolution of a; m >m p for the an- 
nuli and distances from the cluster centre reported from our 
compilation from the literature in Tab.Q] The best match 
between the simulations and the observations is acquired 
for simulations between t — 0At cc and 0.6 t cc , be: we pre- 
dict that the cluster is about half way towards core collapse. 
We therefore conclude that the Arches cluster has not yet 
experienced core collapse but is currently in a pre collapse 
stage. 

In Tab. [1] we have quantified the slope to the low mass 
end of the mass function in 5" to 9" annulus of the Arches 
cluster as consistent w ith Salpeter , wher eas the naive mea- 
surement in Fig. 5 of iKim et all (|2006l ) would results in 
£ m <m p = —3.67 ± 0.14 (with m p — 5Mq), which is unusu- 
ally steep. If this slope would represent the intrinsic Arches 
initial mass function and we adopt a minimum mass of 1 M© 
the observed ~ 2600 stars more massive than ~ 5 A/0 in the 
Arches cluster would result in a total number of more than 
2 x 10 s stars, which is unrealistically high. From an obser- 
vational point of view there are good arguments that the 
low mass end of the mass function is over-estimated, as it 
is plagued by selection effects. One of these effects is the 
artificial correction of missing stars in a crowded field and 
the selection of the three control fields to compensate for 
the background population. In the Keck observations these 
control fields are within about 2.4 pc from the cluster cen- 
tre, which corresponds to ~ 12r CO re- For a King model with 
Wo it 5.2 the control fields would then be located near the 
cluster tidal radius. And since the cluster is about half way 
towards core collapse it is conceivable that the density pro- 
file is described with a King model with Wo ^ 7, in which 
case the control field are part of the cluster halo. 

Due to mass segregation the cluster outskirts will be 
depleted of high mass stars and low mass stars will be over- 
represented (the opposite effect as we discussed for the core 
population) . Correcting the mass function in the cluster core 
with a population taken from near the cluster halo will there- 
fore result in an enormous over correction towards the low 
mass stars, and consequentially result in a steepening of the 
'corrected' mass function. 

One of the control fields (field B of Kim et al, 2006) is 
taken near the location where one expects the tidal tail of 
the cluster in the potential of the Galaxy to pass though. 
The tidal tail is, since it consists of the halo population, also 
likely to be dominated by low mass stars. 

Each of the effects discussed tend to steepen the lower- 
mass end of the mass function, though it is not trivial to 
quantize the effect without a much more detailed study. We 
however, argue that the initial mass function of the Arches 
cluster was probably consistent with Salpeter over the ob- 
served mass range. The observed break in the mass function 
around 5-6 M© and the consequential flattening of the mass 
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function for higher masses is then the result of the dynamical 
evolution of the cluster. 

The initial model which is most comparable to the ob- 
served Arches cluster has a Salpeter initial mass function 
between 1 Mq and 100 Mq and with a reasonably concen- 
trated initial density profile (Wo 4 and Wo 8). 

The break in the mass function in the inner parts of 
the cluster (for r £ 4r core ) appear at m p ~ 2(m), which for 
our simulations is at about 5 Mq . The break in the observed 
mass function in the Arches cluster appears around the same 
mass of m p ~ 5-6 Mq. We performed additional simulations 
with Wo = 9 using a Salpeter mass function down to 0.1 Mq, 
and in this case the break in the core mass function also 
developed around m p ~ 2(m), which for the adopted mass 
function is about 1.0 M0(Figure[5]). Based on these findings 
we argue that the initial mass function in the Arches cluster 
has a lower limit of about 1.0 Mq, as in our simulations that 
reproduce the observations best. 

We predict that other clusters with similar parameter s 
as the Arch es cluster, like Wes terlund 1 dPiatti et al. Il998h . 
NGC 3603 jMoffat et al.ll2004h R13 6 l|Meaburn et alJll982h 
and Quintuplet ( Cotera et al.|[l992T) will show similar char- 
acteristics as Arches. The mass functions in their cores will 
also be rather flat for stars more massive than 5-6 Mq . And 
the mass functions further away from the cluster centre will 
gradually be more like the initial mass function. 
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